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Abstract 

In the context of extended t — J models, with intersite Coulomb interactions of the form 

— n i n ji with n, denoting the electron number operator at site i, nodal liquids are discussed. 

<i,j> 

We use the spin-charge separation ansatz as applied to the nodes of a d-wave superconducting 
gap. Such a situation may be of relevance to the physics of high-temperature superconductivity. 
We point out the possibility of existence of certain points in the parameter space of the model 
characterized by dynamical supersymmetries between the spinon and holon degrees of freedom, 
which are quite different from the symmetries in conventional supersymmetric t — J models. Such 
symmetries pertain to the continuum effective field theory of the nodal liquid, and one's hope is 
that the ancestor lattice model may differ from the continuum theory only by renormalization- 
group irrelevant operators in the infrared. We give plausible arguments that nodal liquids at such 
supersymmetric points are characterized by superconductivity of Kosterlitz-Thouless type. The 
fact that quantum fluctuations around such points can be studied in a controlled way, probably 
makes such systems of special importance for an eventual non-perturbative understanding of the 
complex phase diagram of the associated high-temperature superconducting materials. 



1 Introduction 



The study of strongly correlated electron systems (SCES) is a major enterprise in modern 
condensed matter physics primarily due to high temperature (planar) superconductors, 
fractional Hall conductors, and, more recently, in semiconductor quantum dots. Owing to 
various non-Fermi liquid features of SCES, many believe that the low-energy excitations 
of these systems are influenced by the proximity of a critical Hamiltonian in a generalized 
coupling-constant space. In this scenario, known as spin-charge separation 0, these 
excitations are spinons, holons and gauge fields. 

Important paradigms for SCES are the conventional Hubbard model, or its t — j ex- 
tension, both of which have been conjectured to describe the physics of high-temperature 
superconducting doped antiferromagnets. Numerical simulations of such models , in the 
presence of very-low doping, have provided evidence for electron substructure (spin-charge 
separation) in such systems. 

In ref. ||, an extension of the spin-charge separation ansatz, allowing for a particle- 
hole symmetric formulation away from half-filling, was introduced by writing: 

where the fields z a i obey canonical bosonic commutation relations, and are associated 
with the spin degrees of freedom ('spinons'), whilst the fields ip are Grassmann variables, 
which obey Fermi statistics, and are associated with the electric charge degrees of freedom 
('holons'). There is a hidden non-abelian gauge symmetry SU(2) ® Us(l) in the ansatz, 
which becomes a dynamical symmetry of the pertinent planar Hubbard model, studied in 
ref. f. 

The ansatz (|1|) is different from that of refs. 0, where the holons are represented 
as charged bosons, and the spinons as fermions. That framework, unlike ours, is not 
a convenient starting point for making predictions such as the behaviour of the system 
under the influence of strong external fields. As argued in ||, a strong magnetic field 
induces the opening of a second superconducting gap at the nodes of the <i-wave gap, in 
agreement with recent experimental findings on the behaviour of the thermal conductivity 
of high-temperature cuprates under the influence of strong external magnetic fields [0] . 

In [|J a single-band Hubbard model was used. Such a model should not be regarded 
as merely phenomenological for cuprate superconductors since it can be deduced from 
chemically irealistic multiband models involving both Cu and O orbitals and it has extra 
nearest-neighbour interactions of the form ||: 

2 

Him = -V n i n i U i = Yl C i,i C Ci > ( 2 ) 

<ij> a=l 

as well as longer finite-range hoppings. 

What we shall argue below is that the presence of interactions of the form (§) is crucial 
for the appearence of supersymmetric points in the parameter space of the spin-charge 
separated model. Such points occur for particular doping concentrations. As we shall 



discuss, this supersymmetry is a dynamical symmetry of the spin-charge separation, and 
occurs between the spinon and holon degrees of freedom of the ansatz ([!]). Its appearance 
may indicate the onset of unconventional superconductivity of the Kosterlitz-Thouless 
(KT) type ||, [L0| in the liquid of excitations about the nodes of the d-wave supercon- 



ducting gap ("nodal liquid"), to which we restrict our attention for the purposes of this 
work. 

It should be stressed that the supersymmetry characterizes the continuum relativistic 
effective (gauge) field theory of the nodal liquid. The progenitor lattice model is of course 
not super ■symmetric in general. What one hopes, however, is that at such supersymmetric 
points the universality class of the continuum low-energy theory is the same as that of 
the lattice model, in the sense that the latter differs from the continuum effective theory 
only by renormalization-group irrelevant operators (in the infrared). This remains to be 
checked by detailed studies, which do not constitute the topic of this article. 

In general, supersymmetry provides a much more controlled way for dealing with 
quantum fluctuations about the ground state of a field-theoretic system than a non- 



super symmetric theory [fUJl . In this sense, one hopes that by working in such super- 
symmetric points in the parameter space of the nodal liquid she/he might obtain some 
exact results about the phase structure, which might be useful for a non-perturbative 
understanding of the complex phase diagrams that characterize the physics of the (super- 
conducitng) doped antiferromagnets. As we shall discuss below, to obtain supersymmetric 
points one needs to make specific asumptions about the regime of the parameters of the 
model; from an energetics point of view, such assumptions are retrospectively justified 
by the fact that supersymmetric ground states are characterized by zero energy [[□]], and 
hence are acceptable ground states from this point of view. 

Significant progress towards a non-perturbative understanding of Non-Abelian gauge 
field theories, in four space-time dimensions, based on supersymmetry has been made 
by Seiberg and Witten | I2| . The fact that the spin-charge separation ansatz ([J) of the 



doped antiferromagnet is known to be characterized by such non-Abelian gauge struc- 
ture is an encouraging sign. However, it should be noted that in the case of ref. [0 
extended supersymmetries were necessary for yielding exact results. As we shall discuss 
below, under special conditions for doped antiferromagnets, the supersymmetric points 
are characterized by iV = 1 three-dimensional supersymmetries. Under certain circum- 
stances the supersymmetry may be elevated to N = 2 ||13| , for which it is possible to 



obtain some exact results concerning the phase structure In the present state of the 
udnerstanding of SCES it is a pressing need to have relevant models for which we can 
extract non-trivial exact information. However, for a realistic condensed-matter system 
such as a high-temperature superconductor, even the N = 1 supersymmetry of the super- 
symmetric points is expected to be broken at finite temperatures or under the influence 
of external elctromagnetic fields. Nevertheless, one may hope that by viewing the case 
of broken supersymmetry as a perturbation about the supersymmetric point, valuable 
non-perturbative information may still be obtained. As we shall see, a possible example 
of this may be the above-mentioned KT superconducting properties || that characterize 
such points. 

The structure of the article is as follows: In section 2 we describe briefly the statisti- 



cal model which gives rise to the continuum relativistic effective (2+l)-dimensional field 
theory of the nodal liquid. In section 3 we discuss the properties and (non-abelian gauge) 
symmetries of the spin-charge separation ansatz that characterizes the model. In the next 
section we discuss the intersite Coulomb interactions, which are of crucial importance for 
the existence of supersymmetric points. In section 5 we state the conditions for N=l 
supersymmetry at such points, and describe briefly their importance for yielding super- 
conductivity of Kosterlitz-Thouless type. We conclude in section 6 with some prospects 
for future work. Technical aspects of our work, which may help the non-expert reader to 
follow the arguments presented in the text in a mathematically detailed way, are given in 
two Appendices. 



2 The Model and its Parameters 

In reference || it was argued that BCS-like scenarios for high T c superconductivity based 
on extended t — J models yield reasonable predictions for the critical temperature T c max 
at optimum doping . There it was argued that a pivotal role was played by next-to- 
nearest neighbour and third neighbour hoppings, t' and t" respectively. In particular the 
combination £_ = t' — 2t" determines the shape of the Fermi surface and the nature of 
the saddle points and the associated T™ ax . 

Our aim is to use the extended t — J model studied in []8| in order to discuss the 
appearance of relativistic charge liquids at the nodes of the associated d-wave supercon- 
ducting gap. We will argue that the nodes characterize the model in a certain range of 
parameters. We will demonstrate that at a certain regime of the parameters and dop- 
ing concentration the nodal liquid effective field theory of spin-charge separation exhibits 
supersymmetry. This supersymmetry is dynamical and should not be confused with the 
non-dynamical symmetry under a graded supersymmetry algebra that characterises the 
spectrum of doped antiferromagnets at two special points of the parameter space fT5| . 
We shall also discuss unconventional mechanisms for superconductivity in the nodal liq- 
uid similar to the ones proposed in || PH[ . 

To start with let us describe briefly the extended t — J model used in Ref. ||. The 
Hamiltonian is given by: 

H = P(H hop + Hj + H v )P + PH ll P, (3) 

where: 

Hhop — ~ tij c iot C ja ~~ tij C ia C ja ~ ^ij C ia C ja ) (4) 
(ij) M {ij} 

and (...) denotes nearest neighbour (NN) sites, [...] next-to-nearest neighbour (NNN), 
and {} third nearest neighbour. Here repeated spin (or "colour") indices are summed 
over. The Latin indices i,j denote lattice sites and the Greek indices a = 1,2 are spin 
components. 



(b) 

Hj = J ( Ti^pTj^a — -fliTlj J + J'*^2 Ti,a/3Tj t j3 a , (5) 
(ij) V 4 / m 

2 

with rij = c ta c ia> and Tj Q/ g = c+c^. The quantities J, J' denote the couplings of the 

a=l 

appropriate Heisenberg antiferromagnetic interactions. We shall be interested [|H| in the 
regime where J' « J. 
(c) 

H ^ = llJ2 C ta C ia, (6) 

i 

and \x is the chemical potential, 
(d) 

H v = -V 53 , (7) 

(«> 

This is an effective static NN interaction which, in the bare t — J model, is induced by 
the exchange term, because of the extra magnetic bond in the system when two polarons 
are on neighbouring sites M. Notice that this term, when combined with the Coulomb 
interaction terms in Hj, yields in the effective action a total inter-site Coulomb interaction 
term with coupling 

V to tai = V + 0.25 J (8) 
In ref. [§] the strength of the interaction (|7|) is taken to be: 

V w 0.585 J, (9) 

This is related to the regime of the parameters used in |J, for which the NN hoping 
element satisfies t « J. In fact, for the effective t — j — V model of 0, viewed as an 
appropriate reduction of a single-band Hubbard model, one has the relation: 

At 2 

where U e ff is an effective Hubbard interaction, and J$b is a ferromagnetic exchange 
Heisenberg energy for the single-band model. We have \V'\ ^ \V\ in general, unlike the 
case of the standard Hubbard model with a supplementary intersite Coulomb interaction. 
However, one may consider more general models, in which the above restriction is not 
imposed, and V is viewed as an independent parameter of the effective theory, e.g. 

V J, (11) 

where b is a constant to be determined phenomenologically. Such a situation may arise, for 
instance, in effective models where one considers repulsive on-site Coulomb interactions || 
(e.g. between holes and/or electrons) in addition to the (electron-hole) attractions (0). 
As we shall discuss below, such more general cases turn out to be useful for the existence 
of supersymmetric points in the parameter space of the model. 



(e) The operator P is a projector operator, expressing the absence of double occupancy 
at a site. 

We define the doping parameter < 5 < 1 by 



(m) = i-6, (12) 

d-wave pairing, which seems to have been confirmed experimentally for high-T c cuprates, 
was assumed in ||. A d-wave gap is represented by an order parameter of the form 

A (Jcj = A (cos k x a — cos k y a) , (13) 

where a is the lattice spacing. The relevant Fermi surface is characterised by the following 
four nodes where the gap vanishes: 

±— ,±— J , 14 

2a' 2a J ' V ; 

We now consider the generalized dispersion relation 0, [IB] for the quasiparticles in the 
superconducting state: 



£(£) = y(e (£)-//) +A*(k), (15) 

In the vicinity of the nodes it is reasonable |T^| to assume that \x ~ or equivalently we 
may linearize about /i, i.e. write e (JcJ — /i ~ Vp \q\ where vd is the effective velocity 
at the node and q is the wave- vector with respect to the nodal point. 



3 Non- Abelian spin-charge separation in the t- J model 

As already mentioned in the introduction, it was proposed in ref. || that for the large- 
U limit of the doped Hubbard model the following 'particle-hole' symmetric spin-charge 
separation ansatz occurs at each site i: 



X* = ^ T ,„ 7fti = [ 4 _ C J J f = ( _ V ,t ^ J ( { 22 ^ J ( ("5) 

where the fields obey canonical bosonic commutation relations, and are associated 
with the spin degrees of freedom ('spinons'), whilst the fields ip a ,i, a = 1, 2 have fermionic 
statistics, and are assumed to create holes at the site i with spin index a ('holons'). 
The ansatz (|IBJ) has spin-electric-charge separation, since only the fields ip carry electric 
charge. Generalization to the non-abelian model allows for inter-sublattice hopping of 
holes which is observed experimentally. 

It is worth noticing that the anti-commutation relations for the electron fields c Q ,ct, do 
not quite follow from the ansatz (|16|) . Indeed, assuming the canonical (anti-) commutation 



relations for the z (ip) fields, one obtains from the ansatz ([16]) 



{ c i,i,cl,j} ~ tytM,Aj 

{ C l,h C 2,j} ~ { C 2,i, 4j} ~ 

{ c cM> c a,j} ~ ^« EM^+Wp a = 1,2 no sum over i,j (17) 



13=1,2 



To ensure the canonical anticommutation relations for the c operators we must therefore 
impose at each lattice site the (slave-fermion) constraints 



[ Z i,^h0 + ^P^kil = 1 ( 18 ) 

0=1,2 

Such relations are understood to be satisfied when the holon and spinon operators act 
on physical states. Both of these relations are valid in the large-?/ limit of the Hubbard 
model and encode the non-trivial physics of constraints behind the spin-charge separation 
ansatz (JTBj). They express the constraint of at most one electron or hole per site, which 
characterizes the large-f/ Hubbard models we are considering here. 

There is a local phase (gauge) non-abelian symmetry hidden in the ansatz fll6|) || 
G = SU(2) x Us(l), where SU(2) stems from the spin degrees of freedom, E7g(l) is a 
statistics changing group, which is exclusive to two spatial dimensions and is responsible 
for transforming bosons into fermions and vice versa. As remarked in ||, the Us(l) 
effective interaction is responsble for the equivalence between the slave-fermion ansatz 
(i.e. where the holons are viewed as charged bosons and the spinons as electrically neutral 
fermions [[|) and the slave boson ansatz (i.e. where the holons are viewed as charged 
fermions and the spinons as neutral bosons [T7|, |3|). This is analogous (but not identical) 



to the bosonization approach of |18[] for anyon systems. 

The application of the ansatz (|T6|) to the Hubbard (or t-j models) necessitates a 
'particle-hole' symmetric formulation of the Hamiltonian (|3p, which as shown in 0, is 
expressible in terms of the operators \. In this way, for instance, the NN Heisenberg 
interactions terms become: 

H J = -lT, ^[xixUiXtl (19) 

<ij> 

By making and appropriate Hubbard- Stratonovich transformation on Hj with Hubbard- 
Stratonovich fields A^, we obtain the effective spin-charge separated action for the doped- 
antiferromagnetic model of B: 



H HF = E ( tr [( 8 / J ) A lAi + l^il + CT 3) + \j)4>jVjiUjiip} 



+ 



<ij> 



tr [KziVijUijZj] + h.c.) + ... , (20) 

with the . . . denoting chemical potential terms and NNN hopping terms (the latter are 
essential for the model of M; we shall discuss their effects below). 



This form of the action, describes low-energy excitations about the Fermi surface of 
the theory. The field A y - is matrix valued in 'colour' space; generically it may be expanded 
in components in a canonical basis of 2 x 2 matrices, {1, cr a }, a = 1,2, 3, as follows: 

(A^^^ + i a (a°) ap (21) 

where Greek indices denote 2x2 'colour' indices. 

The quantities Vij and Uij denote lattice link variables associated with elements of 
the £77(2) and Us{l) groups respectively. They are associated || with phases of vacuum 
expectation values of bilinears < 1{Zj > and/or < ip\ ( — t»j (1 + <7 3 ) + A^ipj >. It is 
understood that, by integrating out in a path integral over z and ip variables, fluctuations 
are incorporated, which go beyond a Hartree-Fock treatment. 

The quantity \A X \ is the amplitude of the bilinear < Z{Zj > assumed frozen ||. By an 
appropriate normalization of the respective field variables, one may set |Al| = 1, without 
loss of generality. In this normalization, one may then parametrize the quantity K, which 
is the amplitude of the appropriate fermionic bilinears, as [pj |10|]: 

K= (J|A,| 2 (1-5) 2 ) 1/2 ; l-<5=<£ VW&> , (22) 

a=l 

with 5 the doping concentration in the sample. The quantity |A Z | is considered as an 
arbitrary parameter of our effective theory, of dimensions [energy] 1 ^ 2 , whose magnitude 
is to be fixed by phenomenological or other considerations (see below). To a first ap- 
proximation we assume that A z is doping independent Q. The dependence on J and 5 in 
(f22|) is dictated |TIJ by the correspondence with the conventional antiferromagnetic CP 1 
u-model in the limit 5^0. 

The model of ref. M differs from that of || in the existence of NNN hopping t' and 
triple neighbour hopping t", which were ignored in the analysis of ||. For the purposes 
of this work, which focuses on the low-energy (infrared) properties of the continuum field 
theory of (0), this can be taken into acount by assuming that 

\t i:j \ = t' + = t + 2t+, t + =t' + 2t" (23) 

in the notation of ||. The relation stems from the observation that in the continuum low- 
energy field-theory limit such NNN and triple hopping terms can be Taylor expanded 
(in derivatives). It is the terms linear in derivatives that yield the shift (p3|) of the NN 
neighbour hopping element t. Higher derivatives terms, of the form d x d y are suppressed 
in the low-energy (infrared) limit. 

It is important to note that the model of 0, as well as its extension fl20D , in contrast 
to that discussed in |1, involves only a single lattice structure, with nearest neighbour 
hopping (< ij >) being taken into account, t^. The antiferromagnetic nature is then 
viewed as a property of a 'colour' degree of freedom, expressed via the non-abelian gauge 

1 However, from its definition, as a < ... > of a quantum model with complicated S dependences in 
its couplings, the quantity A z may indeed exhibit a doping dependence. For some consequences of this 
we refer the reader to the discussion in section 6, below. 



structure of the spin-charge separation ansatz flip]). As we shall discuss later, this is very 
important in yielding the correct number of fermionic (holons \P) degrees of freedom in 
the continuum low-energy field theory to match the bosonic degrees of freedom (spinons 
z) at the supersymmetric point. 



4 The Effective Low-Energy Gauge Theory 

4.1 Nambu-Dirac Spinor Representation of nodal Holons 

It is instructive to discuss in some detail the derivation of a conventional lattice gauge 
theory form of the action (|20| ). One first shifts the Ay field: Ay — > Ay- = Ay + ty<73, 
and then assumes that the fluctuations of the Ay field are frozen in such a way that 
only the < A' Q > component is non trivial in the corresponding expansion in terms of 
the Pauli matrices ([H]). This is a variational ansatz that can be justified in the regime 
of the parameters of the statistical model J >> t' + , in which case the dominant Ay- 
configurations (in the path integral) may be taken to be of order J, and thus any effect of 
the (J3 colour structure in the action (|20|) is safely negligible. As we shall discuss in what 
follows, the elimination of the 03 terms from the action fl2"U| ) results in canonical Dirac 
kinetic terms for the fermionic parts of the nodal liquid effective (low-energy) action. 

However, in view of (|23|) , in the model of ||, such an assumption is not valid, given 
that the renormalized hopping parameter, due to NNN and triple neighbour hoppings, 
is of similar order as J. Nevertheless, for our generic purposes in this work we shall 
work in a model where J » t' + . Alternatively, we can assume that the effects of the 
(T3 colour structures can be safely neglected even for the case of the model of ref. ||. 
Such assumptions are retrospectively justified by the fact that the model of || cannot 
yield supersymmetric points even under the above assumption, for other reasons to be 
discussed below. Thus our approach in this paper is to identify the circumstances under 
which deformations of the model presented in || can yield such points in the parameter 
space. 

Notably, the situation J » t' + may be met in the models of Dagotto et al. [||], where 
NNN hopping if is neglected, but where the Coulomb attraction fl7|) is present, in order 
to guarantee the existence of (i-wave superconducting gaps [B5|. Moreover, in the context 



of generalizations of the t — V — J models of Feiner et al. ||, such a situation (c.f. (|10l)) 
is met if one assumes an appropriate attractive V, of opposite sign to the repulsion U e ff, 
but close to it in magntitude (notice that, on account of (23[), in our generalization fo the 
t-V-j model, one should replace t in ( |TU| ) by t' + ). In such a case one has an additional 
large dimensionful scale U e ff, like in the case of the conventional Hubbard model of ||. 

We next remark that in conventional non-abelian gauge theories the fermionic fields 
are usually spinors in the fundamental representation of the gauge group. Let us examine 
under what condition this is feasible in our case. To this end we assemble the fermionic 
degrees of freedom into two 2-component Dirac spinors ||: 

*I,i= (^1 -4) t , *$,,= 4), (24) 



where a in W a i is the colour index. We also consider very weakly coupled SU(2) gauge 
groups, with couplings (?su(2) = Qi « 1- in the weak gauge field approximation, where 
the gauge group element (link) along the fi space-time direction is ~ B®a a + 
0(g 2 ) (with a a , a = 1,2,3 the Pauli matrices), one observes the following mathematical 
identities: 



Tr (ipitp. 



1+J = 



Tr (ViOiV'l+J = *!ri* i+M 
Tr (V^l+J = y\r 3 y i+IM 

Tr (V^l+J = i (-*k 3 ^(ri + iT 2 )* <+/1 + #{1(1 + r 3 )^ +M ) (25) 

where the Pauli matrices r a , a = 1, 2, 3 refer to 'colour' space, and should be distinguished 
from the 03 matrices, which although are 'colour' matrices, they refer to the action (|20|) . 
in which the fermionic degrees of freedom consist of Grassmann variables assembled in 
2x2 matrices. From the last of (|25|), therefore, it becomes evident that the action (|2"0|), 
may be mapped to a conventional lattice action, with spinors (24) in the fundamental 



representation of the 'colour' group, provided that the coupling g 2 << 1 is weak, and in 
addition there is a gauge fixing 0: 



f 3 dx^B 2 = . (26) 

J i 



The weakness of the £77(2) coupling guarantees that a mass gap in the problem is only 
generated by the Us(l) group ||. In the context of the Hubbard model of the coupling 
g 2 of the gauged £77(2) interactions, pertaining to the spin degrees of freedom in the 
problem, is naturally weak, since it is related to the Heisenberg exchange energy J. Given 
that in three space-time dimensions the gauge couplings are dimensionful, with dimensions 
of energy, one may define dimensionless couplings by dividing them with the ultraviolet 
scale of the low-energy theory, which in the model of || is the (strong) Hubbard interaction 
U » J. Thus a dimensionless coupling g 2 ~ J/U << 1 is naturally small in this context. 
A similar situation arises in the context of the effective single-band t-V-j model of ||, 
in the large U e ff >> J limit (c.f. fllOD ). On the other hand, the strong Us(l) coupling 
gi, responsible for mass gap generation for the ho Ions, may be assumed to be of order 
U e ff, since this is the highest energy scale. However, in general for t — j models that 
we consider these relations may not be valid. Still as we shall see below, the ultraviolet 
cut-off of the effective theory, in the regime relevant for supersymmetric points, we are 

2 Note that the requirement for weak gi coupling is essential, given the fact that due to the non-Abelian 
nature of the gauge field, the local gauge fixing B 2 ^ = alone is not sufficient to eliminate dangerous 

terms proportional to a 2 ; this can be easily seen from the Bekker-Hausdorff identity: e 92 "« =1 ' 3 ° Bfl — 



n a =i.3e y Ji f j e 2\92) [v ,0 ]« f4 « (t +...^ wn ^h the commutator being proportional to 02; however such 

terms are of higher order in 52, and hence restriction to weak couplings suffices to yield the conventional 
relativistic gauge form of the effectvie action upon the appropriate gauge fixing. 



interested here, may be up to two orders of magnitude higher than J, thereby allowing 
the Us(l) interactions to be considerably stronger than the SU{2) ones, if one wishes so. 

To generate the conventional Dirac 7-matrix structure for the fermionic action one may 
redefine the spinors in the path integral \I> — > where \I/ are two-component 'coloured' 
spinors, related to the spinors in (p4]) via a Kawamoto-Smit transformation ]20| 



* a (r) = 7o ro • • • 7 2 r2 ^a(r) * a (r) = ¥«(r)( 7 J)*» . . . (ll) ro (27) 

where r is a point on the euclidean lattice, and a = 1, 2 is a 'colour' index, expressing the 
initial antiferromagnetic nature of the system. The generation of Dirac 7 matrices follows 
from identities of the form: 

^A% +fl = yh,AV l+ ,(-ir + ^ + ^ + ^ , A = 1, n, r 3 (28) 

but again the terms proportional to 02 (r 2 ) are problematic and can be eliminated by 
virtue of the gauge fixing (|26|). The 7 matrices appearing in ( |28"D are 2x2 antihermitean 
Dirac matrices on a Euclidean Lattice satisfying the algebra 

{l, , lu} = -2<V (29) 

In terms of the Pauli matrices Oj, i — 1, . . . 3, the 7 matrices are given by 7m = i (03, 01, 02). 
N.B. that fermion bilinears of the form \l/j iQ ,\l/j i( g (i=Lattice index) satisfy 

%, a ^i,P = *i,a*i,/3 (30) 

due to the Clifford algebra (^), and (anti-) hermiticity properties of the 2x2 7 matrices on 
the Euclidean lattice. As we shall see later on, this last identity will be crucial in yielding 
a relativistic form of the effective action for the interacting nodal liquid of excitations in 
generalized Hubbard models. 

We next notice that on a lattice, in the path integral over the fermionic degrees of 
freedom in a quantum theory, the variables \l> and \1/ are viewed as independent. In view 
of this, the spinors ^ in (f23|) may be replaced by as being path integral variables 
on a Euclidean Lattice appropriate for the Hamiltonian system (H). This should be kept 
in mind when discussing the microscopic structure of the theory in terms of the holon 
creation and annihilation operators ip^, ip a , a = 1,2. 

An issue that should be dealt with properly is the appearance of the factors 
(— 1)^+^2 0o+*i)+mi j n (j2^) ) which would prevent the conventional Dirac structure to emerge. 
However, this problem is easily arranged by absorbing such factors in the quantum fluc- 
tuations of the electromagnetic field U(l) em , which are integrated in a path integral 0. 
Notice that by doing so, one does not disturb the form of the bosonic CP 1 parts of 
the effective action, given that the z magnons are electrically neutral. It is understood 

3 A mathematically equivalent, but physically different way, of course, would be to assume a flux 
phase background for the electromagnetic field with flux 7r per lattice plaquette and then consider 
quantum fluctuations about it in a path integral. This would wash out any remnant of the flux phase 
from the effective action, but help in absorbing the above-mentioned factors in a physically irrelevant 
normalization. 



of course, that for the purposes of this work, we shall not be interested further in the 
quantum fluctuations of the electromagnetic field, as their coupling is really much weaker 
than the couplings of the statistical gauge groups under consideration. From now on, the 
electromagnetic interaction will be treated only as an external background. 
The fermionic part of the long-wavelegth lattice lagrangian, then, reads: 

S = ^K'Y^ii-l^Ui,^^ + 

Vi+MU^Vifli] + Bosonic CP 1 parts (31) 

where the Bosonic CP 1 parts denote magnon-field z dependent terms, and are given in 
(p0|). It should be stressed once again that this relativistic form is derived for a weakly- 
coupled SU{2) gauge group, and under a specific gauge fixing. However, in view of the 
gauge invariance, characterizing fl20|) and (|31|), the physical results based on the above 
effective actions, in particular the existence of supersymmetric points in the parameter 
space, of interest to us here, are independent of the the gauge chosen. 

An additional point we would like to make concerns the relativistic form of the action 
(pl|). Although in ( pi]) we did not give explicitly the CP 1 parts, however we have tacitly 



assumed the equality of the effective velocities for spin vs and charge vf (Fermi velocity 
of holes) degrees of freedom. If such an assumption is not made, then the relativistic 
invariant form of of the effective lagrangian is spoiled ||. This can be easily understood 
by the fact that in the effective lagrangian (obtained as a Legendre transform from the 
appropriate Hamiltonian) the (different) velocities vs and vf enter in the derivatives with 
respect to the time variable, e.g. d/vsdt (d/vFdt) in the respective kinetic terms for 
spinons (holons). However, at the supersymmetric points of the nodal liquid, where, 
as we shall discuss later on, the dynamically-generated mass gaps between spinons and 
holons must be equal, the equality vs = vp is essential, otherwise there would be different 
dispersion relations, leading to a difference in mass gaps. These comments should be 
understood in what follows. From now on we shall work in units of the fermi velocity vf- 
The coefficient K' is a constant which stems from the £y— and A^— dependent co- 
efficients in front of the fermion terms in (PU[). An order of magnitude estimate of the 
modulus of (the shifted) • then, which determines the strength of the coefficient K', may 
be provided by its equations of motion. Assuming that the modulus of (the dimensionless) 
fermionic bilinears is of order unity, then, we have as an order of magnitude 

K> ~ (t' + + 1) (32) 



Note that in the regime of the parameters of || t « t + and t + ~ § J for momenta close 
to a node in the Fermi surface, of interest to us here. Thus 



K' ~ 25J/8 (33) 



However, one may even consider more general models, in which K' and the Coulomb 
intersite interaction V are treated as independent phenomenological parameters. 



4.2 Field-Theoretic Treatment of the Constraints 



As discussed in Appendix B, supersymmetrization of CP 1 type models, like the ones 
considered here, requires that the CP 1 constraint be of the form J2a=i \ z a\ 2 — 1- I n our 
case, however, the no-double occupancy constraint, when expressed in terms of the z and 
ty a , a = 1,2, (spinor) fields, with a a 'colour' index, is written as: 

^2[z a z a + ^ a 3 -* a ] = 1 (34) 

a=l 

where (3 = 1/K' 2 , K' is given by (j33|), the 2x2 matrix a% acts in spinor space, and the 
fermions \P are the two- compon ent spinors (0). Equivalently the fermion bilinear terms 
in ([34"D can be expressed in terms of the spinors \& (p7|), which have conventional Dirac 
kinetic terms. This is due to identities of the form (|28|) , extended appropriately to the 
case of coincidence limits in the continuum formalism. The relevant (—1) factors in that 
case may be absorbed in the definition of (3. It is understood that appropriate rescalings 
can be made in the definition of the spinors so as to ensure the canonical kinetic (Dirac) 
term. We have also taken into account that in a euclidean path-integral the variables W 
and \1/ are viewed as independent, which implies that one may redefine W — > \I/ where 
\I/ in later analysis will nevertheless be considered in the conventional way, i.e as ^70- 
Consequently we can interpret the fermion term in the constraint ([33]) as the fermion 
number term. 

The presence of the (non-relativistic) fermion number term in the constraint (|34]) 
appears at first sight to complicate things, since the conventional CP 1 constraint \z\ 2 = 1 
is no longer valid. In fact, as discussed in Appendix B, supersymmetry is compatible with 
the following form of the constraints: 

\z Q \ 2 = 1 , z a m a = (35) 



arising from the superfield version of the CP 1 constraint [13, 2<J. In fact the fermionic 



counterpart of fl35|) can be solved by means of a 'colourless' fermion field X that satisfies 
(on account of the bosonic CP 1 parts of (|35|) ): 

= e a p~zpX , X = e a pz a ^p (36) 

where \1/ Q are the Dirac spinors defined above. To ensure the conventional CP 1 form 
of the bosonic part of the supersymmetric constraints (pB]) from fl34|) we should demand 



(3 « 1, which is satisfied in a regime of the parameters of the theory for which 

K' » K=yfi\& x \{l-8), < 5 < 1 (37) 

For the model of ||, for instance, on account of (p3|), this condition implies that 

\fj/ \A Z \ > 0.32 (1 - 5) , 0<5<1 (38) 

By appropriately rescaling the fermion fields \& to so that in the continuum they 
have a canonical Dirac term, we may effectively constrain the z fields to satisfy the CP 1 
constraint: 



1 2 H ($' — bilinear terms) = 1 

K 



where now the fields \I/ are dimensionful. with dimensions of [energy]. A natural order of 
magnitude of these dimensionful fermion bilinear terms is of the order of K 2 , which plays 
the role of the characteristic scale in the theory, being related directly to the Heisenberg 
exchange energy J. In the limit K' » K fl37|) therefore the fermionic terms in the 
constraint can be ignored, and the constraint assumes the standard CP 1 form involving 



only the z fields ( this being also the case for the model of |9], [L0|, in a specific regime of 
the microscopic parameters). 



As we shall see later, however, the condition (|38 ) alone, although necessary, is not 
sufficient to guarantee the existence of supersymmetric points. Supersymmetry imposes 
additional restrictions, which in fact rule out the existence of supersymmetric points for 
the model of || compatible with superconductivity [J However, this does not prevent 
one from considering more general models in which K' is viewed as a phenomenological 
parameter, not constrained by (|33|). In that case, supersymmetric points may occur for a 
certain regime of the respective parameters. 

However, as a result of the spin-charge separation formalism, there is a different way 
to treat the constraints in a pth integral, which however takes into account the coupling 
of the system to an external elelctromagnetic field, and as such is not apriori relevant 
to the supersymmetric regime. Nevertheless, as we shall discuss in section 6, this will 
be relevant for electric charge transport in the model for which supersymmetry (in the 
absence of external fields) will be argued to play a rather crucial but sabtle role. 

Indeed we observe that the fermion number terms in (|3l| ) may be absorbed in a 
rescaling of the (quantum fluctuations of the) temporal component of the electromagnetic 
field Aq(x, t), which couples (relativistically) only to the spinors \I> (see section 6.2 below). 
Indeed, by implementing the constraint (|3^) in a path-integral via the introduction of a 
Lagrange multiplier field X(x): 

5(\z a \ 2 + /?W a <7 3 * - 1) = / DA(ar)e <A(!B) (*« 2 "+^» < **- 1 ) (39) 

Upon absorbing X(x) in a shift of A (x,t), one obtains from the Maxwell terms in the 
electromagnetic part of the effective action the following combination: 

Hem 3 / o / ox (2di\F i + (di\) 2 ) + standard Maxwell terms (40) 

A(e z /c z ) v ' 

where F 0i is the appropriate components of the Maxwell tensor of the (redefined) elec- 
tromagnetic field, the index i is a spatial index, and repeated indices denote summation, 
The equations of motion for A in the effective action obtained after integrating out, say, 
the z degrees of fredom yield the standard CP 1 model terms BTJ, but also terms of the 



form V^A + 2V t F 0i . One, therefore, may consider a phase in which < X(x) >= const ^ 0, 
provided that the electromagnetic field is chosen as an external one, satisfying Maxwell's 
equations, which is our case. 

The Bosonic part of the constraint, then, implies a mass for the spinons m z oc< 



\{x) > [pl| . The fermionic part on the other hand has the form of a temporal component 



4 We note in passing that in realistic materials superconductivity occurs for doping concentrations 
above 3%, and is destroyed for doping concentrations larger than <5 max ~ 10%. 



of the electric current (see section 6.2 below). The coefficient (3 < X(x) > may be absorbed 
in a shift of the quantum fluctuations of A (x,t). As already stated previously, quantum 
fluctuations of the electromagnetic field will not be of further interest to us here, given 
that we shall treat it only as external background. 

From the above discussion it becomes clear, then, that in either case one maps the 
double occupancy constraint (0) into the standard CP 1 constraint: 

E W 2 = i (4i) 

Q = l 

However, as we have explained above, one cannot really avoid the restriction fl38|), as far as 
the existence of supersymmetric points is concerned, given that any alternative treatment 
would require coupling the system to (supersymmetry-breaking) external electromagnetic 
fields, since otherwise the fermionic parts of fl34|) would be present. As we shall see in 
section 6, though, the alternative treatment of the constraint leads to interesting phases 
of the theory characterized by superconducting elelctric-charge transport. And, then, any 
supersymmetry that might have existed before coupling to elelctromagnetism would play 
an important (but subtle) role in ensuring the existence of superconductivity. 

In addition to the CP 1 constraint, one also encounters the remaining of the constraints 
(Pf), which may also be treated using appropriate Lagrange multiplier fields \2(x), A3 (a;) 
representations for the respective 5-functionals 5(if)\ jV4i)> &{'4 ) i,i' l p2,i)'- 

5{^ 2 )S{^t) ~ J ^rfAs^)^/^^ 2 ^^ 1 ^^ 2 ^^/^^ 3 ^^^^ 2 ^ 1 oc 

/^(^dAa^e^^t" 2 ^^^ 1 ^^^^ 2 ^] (42) 

Above we have expressed the relevant constraint in terms of the spinors fl2~T|), using the 
identity (|28|), appropriately applied to the case of coincident limits in the continuum 
formalism, and absorbed relevant (—1) factors in redefinitions of the lagrange multiplier 
fields Aj, i — 2,3. Notice that the spatial jj, j = 1,2 Dirac matrices are expressed 
in terms of the 2x2 off-diagonal Pauli matrices <jj,j = lm2 as 7^ = iaj. To obtain 
information about the new phases it is necessary to assume < A2(x) >,< As(x) >^ 0. 
We thus observe that the structures in (|42|) resemble terms pertaining to "electric current" 
operators Jj = ^'j^, i = 1, 2 (see sec. 6.2), and as such can be absorbed in the quantum 
fluctuations of the spatial components of the electromagnetic field A(x, t) . 

It should be stressed again that the situation in which the Lagrange multiplier fields 
acquire non-zero vacuum expectation values (vev), < A(x) >, < Aj >^ 0, i = 2,3, cor- 
responds to the selection of a specific ground-state of the system (phase), about which one 
considers quantum fluctuations. There is always the phase in which such vev's are zero, in 
which case one implements the constraints directly on the path-integral correlators, e.g. 
correlation fucntions proportional to ipiip2 are set to zero in this phase. In what follows, 
first we shall resolve the constraints in this latter phase, and later on (section 6) we shall 
discuss the other phases of the model. As we shall later, this phase is characterized by 
spin transport but not electric charge transport, a situation that should be compared 



with the case of the nodal liquids of ref. || in the electrically-neutral-fermion represen- 
tation for spinons. On the contrary, as we shall show in section 6, the phase in which the 
lagrange multiplier vev's are non trivial may yield unconventional superconductivity of 
Kosterlitz-Thouless type || ||]. 

With the above in mind we consider from now on the standard CP 1 constraint involv- 
ing only z fields. By an appropriate normalization of z to z' — ^*_ g the constraint then 
acquires the familiar normalized CP 1 form \z a \ 2 = 1 form. This implies a rescaling of the 
normalization coefficient K in (E0 



K -»■ - = K{1 - 5) ~ A z \ (1 - 5f (43) 

7 

In the naive continuum limit, then, the effective lagrangian of spin and charge degrees of 
freedom describing the low-energy dynamics of the Hubbard (or t — j) model fl20|) of || 
is then: 



C 2 = -Tr 

7 



2 



d l , + ig 2 T a Bl + ig 1 aAz +VD^V (44) 



with z a a complex doublet satisfying the constraint The Trace Tr is over group 



indices, = — ig\a b — ig2T a B a/1 — ^A^, is the gauge potential of the local ('spin') 
SU{2) group, and a M is the potential of the Us(l) group. 

It should be remarked that, we are working in units of the Fermi velocity Vp{— vd) of 
holes, which plays the role of the limiting velocity for the nodal liquid. We stress once again 
that for the nodal liquid at the supersymmetric points we have assumed that vf — vs, 
where vs is the effective velocity of the spin degrees of freedom. The relativistic form of 



the fermionic and bosonic terms of the action ( 44]) is valid only in this regime of velocities. 
This is sufficient for our purposes in this work. Indeed, at the supersymmetric points, 
where we shall restrict our analysis here, the mass gaps for spinons and holons, which may 
be generated dynamically, are equal by virtue of supersymmetry at zero temepratures and 
in the absence of any external fields. Hence it makes sense to assume the equality in the 
propagation velocities for spin and charge degrees of freedom, given that this situation 
is consistent with the respective dispersion relations. This is not true, of course, for 
excitations away from such points. 



5 The NN interaction terms Hy 

We will now discuss the Coulomb-interaction (attractive) terms 

H V = -Vtotal Yl n i n 3 ( 45 ) 
(ij) 

introduced in ref. ||, where V to tai is given in @. With the above discussion in mind 
for the spinors (|24| ) we note that, under the ansatz ([T6D, at a site i the electron number 
operator is expressed, through the Determinant (Det) of the x matrix in (0), in terms 



of the spin, z a , a = 1,2, and charge ip a , a = 1, 2, operators as: 



2 

Deti a/3)i + Det^ a/3ji = ^ (VaV'a + \z a \ 2 ) (46) 



We may express the quantum fluctuations for the Grassmann fields ip a (which now carry 
a 'colour' index a = 1, 2 in contrast to Abelian spin-charge separation models) via: 

where : . . . : denotes normal ordering of quantum operators, and from now on, unless 
explicitly stated, repeated indices are summed over. Since 



V'a.iV'aji) = 1 — ^ , no sum over i 
5 the doping concentration in the sample (0), we may rewrite rtj as 

Tk= (Va| 2 + (l-5)+ : :). 
which in terms of the spinors ^ is given by (c.f. (|2~3|)): 



n< = 2 - <5 + i (*^ 3 *a) . (4S 



where °3 = ( ^1 ) ac ^ S * n ( s P ace "^ me ) s P mor space, and we took into account 



the CP 1 constraint (gT|). 

Consider now the attractive interaction term Hy (45), introduced in ref. ||. We then 
observe than the terms linear in (2 — 5) in the expression for Hy can be absorbed by 
an appropriate shift in the chemical potential, about which we linearize to obtain the 
low-energy theory. We can therefore ignore such terms from now on. 

Next, we make use of the fact, mentioned earlier, that in a lattice path integral the 
spinors W a may be replaced by From the structure of the spinors (f24"D, then, we 
observe that we may rewrite the Hy term effectively as a Thirring vector- vector interaction 
among the spinors \& 



H v = +Y^L £ (* q tA) (*frf*p) (49) 

4 <ij> V J i V / i 



where summation over the repeated indices a,/3(= 1,2), and /x = 0,1,2, with \i = 
a temporal index, is understood. To arrive at fl4T)| ) we have expressed 03 as —270, and 
used the Clifford algerba (|29|) , the off-diagonal nature of the 71^ = io\,i matrices, as well 
as the constraints (ISp. In particular the latter imply that any scalar product between 
Grassmann variables ip a (or r/>t) with different 'colour' indices vanish. 



Taking the continuum limit of (f4~9"l), and ignoring higher derivative terms involving four- 
fermion interactions, which by power counting are irrelevant operators in the infrared, we 
obtain after passing to a Lagrangian formalism 

£v = -j^(v«l^ a ) 2 (50) 

where we have used rescaled spinors, with the canonical Dirac kinetic term with unit 
coefficient, for which the canonical form of the CP 1 constraint (|4l"|) is satisfied. For 
notational convenience we use the same notation \1/ for these spinors as the unsealed 
ones. Although this is called the naive continuum limit, it actually captures correctly the 
leading infrared behaviour of the model. 

We then use a Fierz rearrangement formula for the 7 matrices 

labl^cd = 28ad5bc - Sabred 

where Latin letters indicate spinor indices, and Greek Letters space time indices. The 
Thirring (four-fermion) interactions ( f49[) then become: 

*«7„*«V = -3 (1>A) 2 - 4 £ (¥ a *A*„) (51) 

Notice that this form permits us to use, on account of the identity (j3C|) , either of the 
forms (27) or (^4|) for the spinors or \& in the expression of Hy. It should be noted, 



though, that the canonical Dirac form of the kinetic terms for the spinors is valid only in 
the form (|27|) , which we stick to from now on. 

As mentioned above, in the model of ||, due to the first of the constraints fll8|), the 
mixed colour terms vanish, thereby leaving us with pure Gross- Neveu attractive interaction 
terms of the form: 

£ v = (*A) 2 (52) 

which describe the low-energy dynamics of the interaction (ffij) in the context of the non- 
Abelian spin-charge separation (|16|). It should be stressed that (0) is specific to our 
spin-charge separation model. 

Moreover in the context of the spinors (0), a condensate of the form < > on 

the lattice vanishes because of the constraints (|18D . Such condensates would violate parity 
(reflection) operation on the planar spatial lattice, which on the spinors \l/ is defined to 
act as follows: 

^! (X) -> <7i* 2 (X) , ^ 2 (X) -> <7i*i (X) 

or equivalently, in terms of the (microscopic) ho Ion operxtors ip a , a = 1, 2,: 

1p! (x) Ipl (x) , ^2 (x) -> ~1p\ (x) . 



To capture correctly this fact in the context of our effective continuum Gross-Neveu 
interaction (|52|) the coupling strength must be subcritical, i.e. weaker than the critical 



coupling for mass generation. As discussed in Appendix A, the critical coupling of the 
Gross- Neveu interaction is expressed in terms of a high-energy cut-off scale A as |22fl : 



2 / d 3 q 2g 2 c A 

4 9c / FTT = — — 53 

J Sir^q^ 7i z 

Sa 



where q is a momentum variable and Sa is a sphere of radius A. The divergent g-integral 
is cut-off at a momentum scale A which defines the low-energy theory of interest. For the 
case of interest g 2 = 3 ^ ; t f ' ; on using (|33"1), then, the condition of sub-criticality requires 
that 

A < 77 J . (54) 

which is in agreement with the fact that in all effective models for doped antiferromagnets 
used in the literature the Heisenberg exchange energy J ~ 1000 K serves as an upper 
bound for the energies of the excitations of the effective (continuum) theory. However, 
as mentioned above, to obtain a relativistic gauge theory from the lattice action (|20|) 
one needs the SU(2) interactions to be considerably weaker than the Us(l) interactions, 
responsible for mass generation: the above condition (|54D is also compatible with this, 
provided one identifies the (dimensionful) coupling of the Us(l) interactions with a (high- 
energy) cut-off scale A ~ 77 J. In the context of the effective single-band t-V-j models 
(PH|), for instance, A may be identified with a U e ff » J . 

6 Dynamical Spinon-Holon Symmetry (Supersymme- 
try) in the Nodal Liquid and Potential Phenomeno- 
logical Implications 

6.1 Conditions for N=l Supersymmetry in the nodal liquid 

We turn now to conditions for supersymmetrization of the above continuum theory, i.e. 
conditions for dynamical symmetries between the spinon (boson) and holon (fermion) 
degrees of freedom. Below we shall only outline the main results. Some technical details 
on the formalism are given in |13| and reviewed in Appendix B. Since it has been argued 
that £/s(l) is responsible for dynamical mass generation (and superconductivity) in the 
model of we shall ignore the non-Abelian SU(2) interactions, keeping only the Abelian 
However since the latter argument is not rigorous, it would be desirable to supersym- 
metrise the full group in order to check the phenomenon of dynamical mass generation. 
The extension to supersymmetrizing the full gauge multiplet SU(2) x Us(l) will be the 
topic of a forthcoming work. However we shall still maintain the colour structure in the 
spinors, which is important for the ansatz (|l^) 0. 



As discussed in detail in [13, 23], and reviewed briefly in Appendix B, the conditions 



for N = 1 supersymmetric extensions of a CP 1 o model is that the constraint is of 

5 Ignoring the SU(2) interactions implies, of course, that the 'colour' structure becomes a 'flavour' 
index; however, this is essential for keeping track of the correct degrees of freedom required by supersym- 
metry in the problem at hand |13|, 



the standard CP 1 form fl41|), supplemented by attractive four-fermion interactions of the 
Gross- Neveu type (|5^) , whose coupling is related to the coupling constant of the kinetic 
z-magnon terms of the u-model in a way such as to guarantee the balance between bosonic 
and fermionic degrees of freedom Specifically, in terms of component fields, the pertinent 
lagrangian reads: 

L = gl\D^z a D^z a + i$ + F a F a + 2i(r]^ a z a - W-qiP)] (55) 

where denotes the gauge covariant derivative with respect to the Us(l) field. The 
analysis of [[13], [23| reviewed in Appendix B shows that, upon using the equations of 
motion, 

E^>*) 2 (56) 



F a F 



We thus observe that the N = 1 supersymmetric extension of the CP 1 a model necessi- 
tates the presence of attractive Gross-Neveu type interactions among the Dirac fermions 
of each sublattice, in addition to the gauge interactions. 

In the context of the effective theory (f44"D, (|50|) , discussed in this article, the N = 
1 supersymmetric effective lagrangian (^) is obtained under the following restrictions 
among the coupling constants of the statistical model: 

2 = = 1 ( } 



Note that in the context of the model of ref. ||, for which (p|),(p3|) are valid, the relation 
(|57|) gives the supersymmetric point in the parameter space of the model at the particular 
doping concentration 5 = 5 S : 

(l-5 s ) 2 c^^, 0<5 S <1 (58) 
I 2 1 

According to the discussion in section 4, unbroken supersymmetry (which is valid only 
in the absence of external elelctromagnetic fields) imposes an additional restriction (|38D . 
Then we observe that compatibility of (|58| ) with (|37D , (|38|) requires : 1 — 8 S ^> 1.25, which 
implies that the model of does not have supersymmetric points. 

However, one may consider more general models in which V and K' ~ t' + + J/8 



are treated as independent phenomenological parameters (c.f. (|Tl~|)); in such a case one 
can obtain regions of parameters that characterize the supersymmetric points fl5T|),(|5lf) 
compatible with superconductivity. 

Some comments are now in order: First, it is quite important to remark that in the 
model of 0, where the antiferromagnetic structure of the theory is encoded in a colour 
(non-Abelian) degree of freedom of the spin-charge separated composite electron operator 
(|I|) on a single lattice geometry, there is a matching between the bosonic (z spinon fields) 
and fermionic (\l/ holon fields) physical degrees of freedom, as required by supersymmetry, 
without the need for duplicating them by introducing "unphysical" degrees of freedom |T3|] . 

The gauge multiplet of the CP 1 a- model also needs a supersymmetric partner which 



is a Majorana fermion called the gaugino. As shown in 13 , and reviewed in Appendix B, 



such terms lead to an effective electric-charge violating interactions on the spatial planes, 
given that the Majorana gaugino is a real field, and as such cannot carry electric charge 
(which couples as a phase to a Dirac field). These terms can be interpreted as the removal 
or addition of electrons due to interlayer hopping. 

Indeed, the gaugino r\ terms in the supersymmetric lagrangian fl55|) have the form: 

and hence may be viewed heuristically as constituting a Majorana-spinor representation 
of the absence of spin and charge at a site of the planar lattice system To understand this, 
the reader is advised to make a comparison with the Grassmann XiX ■> representation of 
a Wilson line ('missing spin' S ) in the treatment of static holes in refs. |25, ||: 

d X ] d X e- lS 'J ■*E 4 (-i) < xt»ao(i,t) ( 60 ) 

where a is the temporal component of the gauge potential of the CP 1 a- model, describing 
spin excitations in the antiferromagnet. From this point of view, the existence of N — 1 
supersymmetry in the doped antiferromagnets necessitates interplanar couplings, through 
hopping of spin and charge degrees of freedom (electrons) across the planes. In view of 
( [571) such interlayer hopping is suppressed by terms of order \fj . 

Another important point we wish to make concerns the four-fermion attractive Gross- 
Neveu interactions in fl55D ,(|56f). As discussed in detail in f24|, if the coupling of such 



terms is supercritical, then a parity-violating fermion (holon) mass would be generated in 
the model. However, the condition (|54D , which is valid in the statistical model of interest 
to us here, implies that the respective coupling is always subcritical, and thus there is no 
parity-violating dynamical mass gap for the holons, induced by the contact Gross-Neveu 
interactions. This leaves one with the possibility of parity conserving dynamical mass 
generation, due to the statistical gauge interactions in the model || |24j. 

A detailed analysis of such phenomena in the context of our CP 1 model is left for 
future work. For the present, however, we note that in N = 1 supersymmetric gauge 



models, supersymmetry-preserving dynamical mass is possible 0, In fact, as 

discussed in [ETJ, although by supersymmetry the potential is zero, and thus there would 
naively seem that there is no obvious way of selecting the non-zero mass ground state over 
the zero mass one, however there appear to be instabilities in the quantum effective action 
in the massless phase, which manifest themsleves through instabilities of the pertinent 
running coupling. The opening of such a fermion mass gap has been associated with the 
existence of a non-trivial infrared fixed point of the renormalizat ion-group flow, which 
implies non-fermi liquid behaviour [ES . 



From a physical point of view, such a phenomenon would imply that, for sufficiently 
strong gauge couplings, the zero temperature liquid of excitations at the nodes of a d wave 
superconducting gap would be characterized by the dynamical opening of mass gaps for the 
holons. At zero temperature, and for the specific doping concentrations corresponding to 
the supersymmetric points, as advocated above, the nodal gaps between spinon and holons 
would be equal, in agreement with the assumed equality of the respective propagation 



velocities vp = vs, which yielded the relativistic form of the effective continuum action 
([44]) of the nodal excitations at the supersymmetric points. 

The opening of a nodal mass gap, due to the Us(l) gauge interactions, would imply 
a breaking of the fermion number (global U(l)) symmetry, and thus superconductivity 
upon coupling the system to external electromagnetic fields, according to the scenario of 
H |3|, which is reviewed briefly below for the benefit of the non-expert reader. 

6.2 Kosterlitz-Thouless Realization of Superconductivity in the 

577(2) <8> U s (l) model 

This section is mainly a review of results that appear in the literature regarding the 
model 0, ||, [24H - It mainly serves as a comprehensive account of the various delicate 



issues involved, which play a very crucial role in the underlying physics. It is primarily 
addressed to the non-experts in the area. Only the basic results will be presented; the 
interested reader may then find the relevant details in the published literature. 

An important issue in the effective gauge theory SU(2) <g> Us{l) model is the existence 
of a global conserved symmetry, namely the fermion number, which is due to the electric 
charge of the fermions The corresponding current is given by 

2 

J„ M =0,l,2. (61) 

a=l 

This current generates a global Ue(1) symmetry, which after coupling with external elec- 
tromagnetic fields is gauged. In this sense the holon current (|6~1~D coincides with the charge 
transport properties of the system. 

Some discussion is in order at this point. The association of the current J M ([H]) with 
an electric current for holons comes about due to the similarity of the form of the spinors 
([24]) with the conventional Nambu spinors appearing in the BCS Hamiltonian for super- 
conductivity. Indeed, for the benefit of the reader we remind that in such a case the 
electron operators c a are assemblied, in a particle-hole formalism, into two component 
spinors (of,cj), and the resulting Hamiltonian couples in a gauge invariant way to an 

external electromagnetic potential A by making the standard substitution of the momen- 
tum operator p — > p — -A. The only difference in our nodal liquid case is that the holon 
spinors (^) come in two 'colours' and, as contrasted to the generic BCS case, the prob- 
lem is relativistic due to the restriction in the nodal excitations. Thus, at the level of 
the continuum effective action of the nodal excitations, the coupling to electromagnetic 
potentials is straightforward by extending the (statistical) gauge covariant derivatives in 
the Dirac kinetic terms ( [44]) to incorporate the electromagnetic potential coupling terms 

J d 3 x-j2^ a j,A^ a (62) 



where c is the light velocity and e is the absolute value of the electron charge (for holon 
excitations the charge is +e, for electron — e; in our problem here we concentrate in the 



holon current). The resulting nodal holon electric current is given by differentiation with 



respect to A M , i.e. by the expression (plf) . 

Before discussing superconducting properties of the system we should remark that, as 
a result of the constraints (|18]) and the non-diagonal nature of the 7,-, z = 1,2 matrices, 
the spatial components of the current (^T|) vanish, but the temporal component (charge 
density) is non trivial. Moreover, given that the constraints ([18]) do not concern the 
spinons z, this means that there is a phase of the nodal liquid in which there is no charge 
transport, but only spin transport. The non-trivial 'spin current' may be thought of as 
given by J^ pin ~ zd^z. This situation should be compared with the corresponding phase 
in nodal liquids in the approach of ref. [j5|, where the spinons are represented as electrically 
neutral fermions. 

However, in our model there are other possibilities, leading to more complicated 
phases, as we shall discuss now. These possibilities are realized by implementing the con- 
straints flrgp via appropriate lagrange multipliers in the path integral over the fermionic 
variables ^ , ip, as we discussed in section 4 (c.f. (p9|),(^2|)). Expressing the products iftifa 
(and their conjugates) as spatial components of the current (|6i|) , then, one may assume 
a specific ground state in which the appropriate lagrange multipliers for the constraint 
V'i'02 ~ (and hermitean conjugate) acquire non-zero vacuum expectation values that 
may be absorbed by appropriate shifts of the corresponding spatial components of the 
electromagnetic potential A(x,t) coupled to the current J. As we have already discussed 
in section 4, a non-trivial vacuum expectation value for the lagrange multiplier X(x) of the 
last of the constraint (|T8| ) will yield mass terms for the z magnons, whilst the fermionic 
part of the constraint may be absorbed by an appropriate shift of the temporal component 
of the electromagnetic potential. This procedure breaks supersymmetry explicitly but, 
as we shall argue now, the existence of supersymmetry before coupling to external elec- 
tromagnetism is crucial in implying superconducting properties after coupling to external 
fields. 

In this framework, the constraints ([IB] ) no longer apply in the path integral, and non- 
vanishing spatial compontents of the electric current, J, appear. It should be remarked 
that in such a case the mixed colour terms in ([51]) do not vanish, and hence the resulting 
effective lagrangian breaks supersymmetry explicitly. This was to be expected, anyhow, 
from the the very presence of external (non supersymmetric) electromagnetic fields. How- 
ever, given that the coupling of such contact four fermion interactions is subcritical (c.f. 
(^),(Q)), such interactions are irrelevant operators in a renormalization-group sense, and 
hence the universality class of the theory (in the infrared) can still be determined using 
the supersymmetric version of the theory in the absence of any external fields (which also 
satisfies the additional restriction (p8|)). As we shall argue below, this more general phase 
is important in that it yields unconventional superconductivity for the nodal liquid. 

To this end, we remark that in the absence of external electromagnetic potentials, the 
symmetry Ue{1) is broken spontaneously in the massive phase for the fermions This 
can be readily seen by considering the following matrix element (see figure [I]): 

S a =<B°\J„\0>, a = 1,2,3 ; J„ = # 7 „* (63) 



Figure 1: Anomalous one-loop Feynman matrix element, leading to a Kosterlitz-Thouless- 
like breaking of the electromagnetic U em (l) symmetry, and thus superconductivity, once 
a fermion mass gap opens up. The wavy line represents the SU(2) gauge boson B 3 , 
which remains massless, while the blob denotes an insertion of the fermion-number current 
J M = \I / 7 M \1/. Continuous lines represent fermions. 



As a result of the colour group structure only the massless gauge boson of the SU(2) 
group, corresponding to the 03 generator in two-component notation, contributes to the 
graph. The result is |9|, [29| : 



S =< Bl\J u \0 >= ( S gnM) W -^|: (64) 



where M is the parity-conserving fermion mass (or the holon condensate in the context 
of the doped antiferromagnet). In our case this mass is generated dynamically by means 
of the Us(l) interactions, as we discussed above, provided its coupling constant is suffi- 
ciently strong. The result (|64]) is exact in perturbation theory, in the sense that the only 



modifications coming from higher loops would be a multpilicative factor fZnTp) 011 ^ ne 
right hand side, with Tl(p) the -B^-gauge-boson vacuum polarisation function |29[ 



As discussed in || g9[, the B colour component plays the role of the Goldstone boson of 
the spontaneously broken fermion-number symmetry. If this symmetry is exact, then the 
gauge boson B^ remains massless. This is crucial for the superconducting properties |§, 
given that this leads to the appearance of a massless pole in the electric-current two-point 
correlators, the relevant graph being depicted in figure §. This is the standard Landau 
criterion for superconductivity. 

It can be shown || that in the massive- fermion (broken 577(2)) phase, the effective 
low-energy theory obtained after integrating out the massive fermionic degrees of freedom 
assumes the standard London action for superconductivity, the massless excitation being 
defined to be the dual of B^: 

= e^ p d u B 3 p (65) 

All the standard properties of superconductivity, Meissner effect (strongly type II ||), 
flux quantization and infinite conductivity, follow then in a standard way after coupling 
to external elelctromagnetic potentials, provided the excitation cf) (and, hence, B 3 ) is 
exactly massless. 




Figure 2: The lowest-order contribution to the electric current- current correlator 
(0| Jfj,{p) J v {— p) 1 0) . The blob in the propagator for the gauge boson B^ indicates fermion 
loop (resummed) corrections. The blob in each fermion loop indicates an insertion of the 



current J^. 

However, it is known |], [TU], |3| that superconductivity is of a Kosterlitz-Thouless (KT) 
type superconductivity, not characterized by a local order parameter. Let us briefly review 
the arguments leading to this ||. The neutral parity-invariant condensate < ^/i^i — 
^2^2 >, generated by the strong Us{l) interaction, is invariant under the U(l) <S> Ue(1), 
as a result of the t% coupling of B^ in the action, and hence does not constitute an order 
parameter. This is a characteristic feature of our gauge interactions. Putative charge 2e or 
— 2e order parameters, like the pairing interactions among opposite spins in the statistical 
model of |], [J, e.g. < >, < ^i 1 !^ > Q will vanish at any finite temperature, in the 

sense that strong phase fluctuations will destroy the vacuum expectation values of the 
respective operators, due to the Mermin- Wagner theorem. Even at zero temperatures, 
however, such vevs yield zero result to any order in perturbation theory trivially, due to 
the fact that in the context of the effective B^ gauge theory of the broken SU(2) phase, 
the gauge interactions preserve 'flavour'. For a more detailed discussion on the symmetry 
breaking patterns of (2 + l)-dimensional gauge theories, and the proper definition of order 
parameter fields, we refer the reader to the literature [J. Thus, from the above analysis 
it becomes clear that gap formation, pairing and superconductivity can occur in the above 
model without implying any phase coherence. 

6.3 Instantons and the fate of Superconductivity in the 577(2) 
U s (l) model 

An important feature of the non-Abelian model is that, due to the non-Abelian symmetry 
breaking pattern SU(2) — ► U(l), the abelian subgroup U(l) G SU(2), generated by the 
cr 3 Pauli generator of SU(2), is compact, and may contain instantons [Q, which in three 
space-time dimensions are like monopoles, and are known to be responsible for giving a 

6 In four-cmponent notation, such fermionic bilinears correspond to < '5 75 '5 >, < 'I' 75 '5 >, considered 
in §. 



small but non-zero mass to the gauge boson B z 



m B 3 ~ e 2 So (66) 

where Sq is the one-instanton action, in a dilute gas approximation. Its dependence on 
the coupling constant g% = gsu{2) is well known |30| : 

C COIlst (a7\ 

For weak coupling g<i the induced gauge-boson mass can be very small. However, even 
such a small mass is sufficient to destroy superconductivity, since in that case there is no 
massless pole in the electric current-current correlator. In a breakdown of supercon- 



ductivity due to instanton effects has been interpreted as implying a "pseudogap" phase: 
a phase in which there is dynamical generation of a mass gap for the nodal holons, which, 
however, is not characterized by superconducting properties. 

The presence of massless fermions, with zero modes around the instanton configu- 



ration, is known |30| to suppress the instanton effects on the mass of the photon, and 
under certain circumstances, to be specified below, the Abelian-gauge boson may remain 
exactly massless even in the presence of non-perturbative effects, thus leading to super- 
conductivity, in the context of our model. This may happen |30| if there are extra global 
symmetries in the theory, whose currents connect the vacuum to the one -gauge-boson 
state, and thus they break spontaneously. This is precisely the case of the fermion num- 
ber symmetry considered above |3l], |2"5fl . In such a case, the massless gauge boson is the 
Goldstone boson of the (non-perturbatively) spontaneously broken symmetry. However, 
in our SU(2) £g> Us(l) model [^, as a result of the (strong) Us(l) interaction, a mass 
for the fermions is generated, and hence there is no issue of fermion zero modes in this 



case. The analysis of the low energy effective theory presented in ||, |24| is based on a 
Wilsonian treatment, where massive degrees of freedom are integrated out in the path 
integral. This includes the gapful fermions and the massive SU(2) gauge bosons. The 
resulting effective theory, then, which encodes the dynamics of the gapped phase, is a 
pure gauge theory U{1) G SU(2), and the instanton contributions to the mass of are 
present, given by in the one-instanton case. Thus, it seems that, generically, in the 
context of the SU{2) ® Us(l) of ref. 0, the nodal gap is actually a pseudogap. 



6.4 Instantons and Supersymmetry 

We now remark that Supersymmetry is known |JIJ to suppress instanton contributions. 
For instance, in certain N = 1 supersymmetric models with massless fermions, considered 
in ref. the instanton-induced mass of the Abelian gauge boson is given by: 

gauge boson ~ 6 (68) 

which is suppressed compared to the non-supersymmetric case (|B6j). 

N = 2 supersymmetric theories in three space-time dimensions constitute additional 
examples of theories where the abelian gauge boson remains exactly massless, in the 



presence of instantons |30], [!![]. Such theories have complex representation for fermions, 
and hence are characterized by extra global symmetries (like fermion number). In view of 
our discussion above, such models will then lead to Kosterlitz-Thouless superconductivity 
upon gauging the fermion number symmetry. 
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Figure 3: A possible scenario for the temperature-doping phase diagram of a charged, rela- 
tivistic, nodal liquid in the context of spin-charge separation. At certain doping concentra- 
tions (Sss) there are dynamical supersymmetries among the spinon and holon degrees of 
freedom, responsible for yielding thin "stripes" in the phase diagram (shaded region) char- 
acterized by Kosterlitz-Thouless (KT) superconductivity without a local order parameter. 
The diagram is conjectural at present. It pertains strictly to the nodal liquid excitations 
about the d-wave nodes of a superconducting gap, and hence, should not be confused with 
the phase diagram of the entire (high-temperature) superconductor. 



In this respect, the supersymmetric points (p7D,(p8[) for which such instanton effects 
are argued [E3 to be strongly suppressed in favour of KT superconductivity, as reviewed 



above, would constitute "superconducting stripes" in the temperature-doping phase di- 
agram of the nodal liquid (see fig. ||]) 0. Theoretically, the stripes should have zero 
thickness, given that they occur for specific doping concentrations fl5?|) , (p8|) . However, in 
practice, there may be uncertainties (due to doping dependences) in the precise value for 
the parameter A z entering (|57D, (|38|) which might be responsible for giving the supercon- 
ducting stripe a certain (small) thickness. A detailed analysis of such important issues is 
still pending. It is hoped that due to supersymmetry one should be able to discuss some 
exact analytic results at least for zero temperatures. 

We also remark that in supersymmetric theories of the type considered here and in 



ref. [13], it is known [30 that supersymmetry cannot be broken, due to the fact that 



7 It should be stressed that the term "stripe" here is meant to denote a certain region of the 
temperature-doping phase diagram of the nodal liquid and should not be confused with the stripe struc- 
tures in real space which characterizes the cuprates at special doping concentrations. 



the Witten index (— 1) F , where F is the fermion number, is always non zero. Thus, 
in supersymmetric theories the presence of instantons should give a small mass, if at 
all, in both the gauge boson and the associated gaugino, However, in three dimensional 
supersymmetric gauge theories it is possible that supesymmetry is broken by having the 
system in a 'false' vacuum, where the gauge boson remains massless, even in the presence 
of non perturbative configurations, while the gaugino acquires a small mass, through non 
perturbative effects. The life time, however, of this false vacuum is very long and 
hence superconductivity can occur, in the sense that the system will remain in that false 
vacuum for a very long period of time, longer than any other time scale in the problem. 



6.5 Some Comments on Supersymmetry Breaking at finite tem- 
peratures 

So far, our discussion was restricted to zero temperature. At any finite temperature, 
no matter how small, supersymmetry is explicitly broken, and thus the supersymmetric 
points should be viewed as quantum critical points. However, the breaking of supersymme- 
try is associated with different boundary conditions between fermionic and bosonic degrees 
of freedom, and, although the vacuum energy is no longer zero, however a detailed analy- 
sis should be made in order to determine whether the equality of mass gaps between the 
nodal spinons and holons at the supersymmetric points is lifted by temperature-dependent 
corrections. In the context of a supersymmetric theory this issue can be tackled by means 
of "thermal superspace" methods, which have been developed recently in the context of 



particle-physics models [32]. The generic result of such analyses seems to be that the 
mass degeneracy among the superpartners is lifted at the level of the mass of the various 
thermal modes, the corresponding lifting being proportional to the temperature. The 
thermal superspace method can be applied to the present model as well, however this 
falls beyond the scope of the present article and is thereby left for a future work. 

Moreover, as the crude analysis of indicates, the nodal gaps would disappear at 
temperatures which are much lower than the critical temperature of the (bulk) <i-wave 
superconducting gap. For instance, for a typical set of the parameters of the t — j model 
used in ||, the nodal critical temperature is of order of a few mK, which is much smaller 
than the 100 K bulk critical temperature of the high temperature superconductors. The 
application of an external magnetic field in the strongly type II high-temperautre super- 
conducting oxides, which is another source for explicit breaking of the potential super- 
symmetry, enhances the critical temperature || up to 30 K, thereby providing a potential 
explanation for the recent findings of 0, according to which plateaux in the thermal 
conductivity as a function of the external magnetic field indicate the opening of a gap at 
the o?-wave nodes. 

We now remark that, if such situations with broken supersymmetry are viewed as cases 
of perturbed supersymmetric points, then one might hope of obtaining non-perturbative 
information on the phase structure of the liquid of nodal excitations in spin-charge sepa- 
rating scenaria of (gauge) high-temperature superconductors. This may also prove useful 
for a complete physical understanidng of the entire phenomenon, including excitations 
away from the nodes. 



7 Conclusions 



From the above discussion it is clear that supersymmetry can be achieved in the effec- 
tive continuum field theories of doped antiferromagnetic systems exhibiting spin-charge 
separation only for particular doping concentrations (cf. ( |58|) ,(|38|)). One's hope is that 
the ancestor lattice model will lie in the same universality class (in the infrared) as the 
continuum model, in the sense that it differs from it only by the action of renormalization- 
group irrelevant operators. This remains to be checked by explicit lattice calculations. 
We should note at this stage that this is a very difficult problem; in the context of four- 
dimensional particle-physics models it is still unresolved [Q. However, in view of the 
apparent simpler form of the three-dimensional lattice models at hand, one may hope 
that these models are easier to handle. 

By varying the doping concentration in the sample, one goes away from the super- 
symmetric point and breaks supersymmetry explicitly at zero temperatures. At finite 
temperatures, or under the influence of external electromagnetic fields at the nodes of the 
d-wave gap, supersymmetry will also be broken explicitly. Therefore, realistic systems 
observed in nature will be characterized by explicitly broken supersymmetries even close 
to zero temperatures. However there is value in deriving such supersymmetric results in 
that at such points in the parameter space of the condensed-matter system it is possible 
to obtain analytically some exact results on the phase structure of the theory. Supersym- 
metry may allow for a study of the quantum fluctuations about some exact ground states 
of the spin-charge separated systems in a controlled way. Then one may consider perturb- 
ing around such exact solutions to get useful information about the non-supersymmetric 
models. 

We have argued that such special points will yield new phases for the liquid of exci- 
tations about nodal points of the d-wave superconducting gaps, which include a phase in 
which there is only spin transport but not electric current transport, as well as a phase 
in which there are Kosterlitz-Thouless type superconducting "islands" in a temperature 
doping phase diagram of the nodal liquid, upon the dynamical generation of holon-spinon 
mass gaps (of equal size). The latter property is due to special properties of the super- 
symmetry, associated with the suppression of non-perturbative effects of the (compact) 
gauge fields entering the spin-charge separation ansatz (p]). This, of course, needs to be 
checked explicitly by carrying out the appropriate instanton calculations in the spirit of 



the non-perturbative modern framework of 12|. At present, such non-perturbative effects 



can only be checked explicitly in three dimensions for highly extended supersymmetric 
models pl| . It is, however, possible that some exact results could be obtained at least 
for the N = 2 supersymmetric models which may have some relevance for the effective 
theory of the nodal liquid at the supersymmetric points JTBJ. Then, one may get some 



useful information for the N = 1 models studied here by viewing them as supersymmetry- 
breaking perturbations of the N = 2 models. Such issues remain for future investigations, 
but we hope that the speculations made in the present work provide sufficient motivation 
to carry out research along these directions. 
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Appendix A 



Renormalization aspects of Four-Fermi Theories in fewer than four space- 
time dimensions 



In this appendix we shall review briefly the renormalization-group approach to rela- 
tivistic theories with four-fermion interactions in fewer than four space-time dimensions. 
Below we shall outline only the basic results. For further details we refer the interested 



reader to ref. p2 |. 

We shall use as our pilot theory a three-dimensional model, with four-component 
spinors, containing Gross-Neveu interactions. The lagrangian is given by: 



9 

2N 



(69) 



where i = 1,2, ... N, N is a fermion species ('flavour') number, which is assumed large, 
and repeated indices % denote summation. 

Linearizing, by means of a Hubbard- Stratonovich scalar (a) field the four-fermion 
interactions yields: 



C 



1 

27 2 



1 



(70) 



By naive power counting the four-fermion terms are irrelevant non-renormalizable 



operators. However, the basic observation [^] was that in the large N limit the ultraviolet 
behaviour of the fermion propagator is softened in such a way that the scaling dimension 
of the composite operator ifjip changes from its naive dimension, so that the four-fermion 
interactions become renormalizable. 

This can be seen as follows: from ( [70|) we observe that the tree level scalar propagator 
is given by 



G<°>(p) 



9 



(71) 



Consider now the one-loop fermion-vacuum polarization graph. Assume that a fermion 
mass m is generated dynamically, m can be determined self-consistently by a Schwinger- 
Dyson approach ||22|| . For our purposes the details of the derivation will be omitted. The 



result for the one- loop vacuum polarization graph is [22 



n( P ) 



-Tr 



d 3 k 



1 
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-2F(p) 



(72) 



In the large iV limit the loop graphs can be resummed. To leading order in 1/N 
expansion the dressed propagators can be expresssed as 



G(p) 



g~ 2 + 2F(p) 



(73) 



From this epxression it is obvious that in the ultraviolat limit p — > A >> m, where 
A an ultraviolet (high-energy) cut-off, the behaviour of the propagator G(p) is such that 
G(p) ~ i, which implies that the scaling (mass) dimension of the field a in the ultraviolet 
regime is [o] uv = 1. From the action (|70|) it is obvious that the field a is equivalent to 
the composite field ipip, as far as scaling dimension is concerned. This implies then that 
the mass dimension of the four-fermion Gross-Neveu operator is 

= 2 . (74) 

uv 

This guarantees the renormalizability of the theory, since the pertinent operators be- 
came relevant in a renormalizat ion-group sense: if one computes the effective quantum 
corrections to the four-fermion scattering amplitude in the resummed -k approximation, 
then the result for the renormalized coupling is given by G(p) in (|73|) which for p » m 
scales like 2 

9l = (75) 

showing that the effective interaction grows strong for high momenta for real g (attractive 
four fermion interactions in our notation). Notice that for the four-fermion theory the 
renormalizability concerns the ultraviolet (high-energy) regime. There is an UV stable 
non-trivial fixed point in the theory and the associated critical exponents can be computed 
within the -k expansion up to order [ 22 1 . Such computations have also been compared 



successfully with corresponding results from lattice simulations. 

A Schwinger-Dyson analysis for mass generation |22[] leads, in the large A" limit, to 
the following gap equation: 

tm = A92 lwf q \ q ^ +m -) (76) 

2 2 

with m the dynamically generated mass, and t = 9 ~2 C , and the critical coupling g c is 
defined through 1 = Ag 2 c J §4,. 

Thus, mass generation occurs only for positive four-fermion couplings (attractive) 
Gross-Neveu interactions, which are stronger than a given critical value. In a renormal- 
ization group sense the repulsive interactions are irrelevant, becoming weaker and weaker 
as one lowers the momenta. 

Since the supersymmetric version of the CP a-model, of interest to us here, contains 
- as we discuss in Appendix B - both Gross-Neveu and Thirring [tp^^ip) 2 interactions in 
it component form ]13| , we turn next our attention to a brief review of a renormalization 



group study of such mixed models. 

Such models have been discussed in the literature E^|, with the conclusion that it is 
mainly the Gross-Neveu interactions which determine the critical behaviour, in a large A" 
framework. Let us review the situation briefly. The lagrangian is given by: 

2 7,2 

C = ^ fhl> + (m) 2 + — (^ 7 ^) 2 (77) 



where i — 1, 2, . . . N, N is a fermion species ('flavour') number, which is assumed large, 
h 2 is the coupling of the Thirring interactions, and, as before, repeated indices denote 
sunmmation. 

The Thirring interactions become renormalizable in the UV, just as the Gross Neveu 
ones, which can be proven in a similar way to the Gross Neveu interactions above, i.e. 
by linearizing the Thirring interaction by a Habbard-Stratonovich vector field A^. The 
vector interactions are viewed as gauge fixed interactions with a bare propagator fl22" | 

A^(p) = h 2 [5^ - Pf,p u /p 2 ) + gauge - fixing terms (78) 

The dressed (in 1/N expansion ) vector propagator is modified 

A Mi ,(p) = h -2 + F (pj (V - P^Pu/p 2 ) + gauge - fixing terms (79) 

where F{p) has been defined in (|72|). 

In the ultraviolet regime p — > A the scaling mass dimension 

of the vector field ~ ip^^ip is again one, leading to a renormalizable Thirring 
interaction in the ultraviolet. 

A detailed analysis |22[ of the critical behaviour in this combined Gross-Neveu and 
Thirring model shows that the critical behaviour is driven by the UV fixed point of the 
Gross-Neveu terms. Moreover, repulsive Thirring terms cannot lead to dynamical mass 
generation, and thus do not affect the critical (fixed point) behaviour of the theory. 

This analysis implies that, up to irrelevant operators in a renormalization group sense, 
from the various four-fermion contact interactions in our effective theory, the Gross-Neveu 
type interactions appearing in fl52|) are the only ones that could affect the universality 
class of the model, leading to a non trivial Ultraviolet stable fixed point. However, in 
the context of the planar condensed matter systems with relativistic fermions we are 
discussing here, the Gross-Neveu four-fermion contact interactions are sub-critical (c.f. 
([54])), and hence irrelevant operators in the infrared (low energy) limit. In our systems 
it is the gauge-field-holon interactions that grow strong for low momenta and are thus 
relevant in a renormalization group sense. This point has been discussed in detail in 



where we refer the interested reader. In fact, such interactions have been argued to be 
responsible for a non-fermi liquid behaviour of the pertinent relativistic liquids. 

Appendix B 

N=l Supersymmetric CP 1 a-models in (2+l)-dimensions 

In this Appendix we shall be interested in discussing briefly the formalsim underlying 
supersymmetrization of a CP 1 model coupled to Dirac fermions: 



£2 = gl\{d^ - a^)z\ 2 + MDtffjS! 



(80) 



where now = d^—ia^, g\ has dimensions of mass, a M is the Us(l) ('fractional statistics') 
field. For simplicity we consider as a gauge interaction that of a standard Us(l) Abelian 
gauge theory. 

We consider the standard CP 1 constraint: 

E = 1- (81) 

a=l 

As we shall discuss immediately below, this form of the constraint can be supersym- 
metrized. 

We now proceed to the supersymmetrization of the model ( 80|) with the constraint 
(pl|). Below we shall outline only the main results. For details we refer the reader to ref. 



fl3| , p3(1 and references therein. The main idea behind such a supersymmetrization is to 
view the magnons z as supersymmetric partners of the holons 

The basic "matter" multiplet of N=l supersymmetry in three space-time dimensions, 
can be written in terms of a scalar superfield as 

$ = (/> + X + (1/2)0GF (82) 

which contains a real scalar field, 0, a Majorana spinor x and a real auxiliary field F. We 
consider complex superfields 

Z = (l/v / 2)($i + z$ 2 ) = z + m + (1/2)99F (83) 

which contain a complex scalar, z = (l/v / 2)(0i + ^2)5 a Dirac spinor, \1> = (l/v / 2)(xi + 
1X2)1 and a complex auxiliary field, F = (l/v / 2)(i r i + iF^). The supersymmetry transfor- 
mations read, 

5 s z = ^ 

8s* = -i'f£d ll z + £F 

5 S F = -i£ jd^! (84) 

and the supersymmetric invariant lagrangian is given by the highest component (99) of 
the superfield DZ*DZ, where 

D a = JL-i(P9) a (85) 

is the supersymmetry covariant derivative. 

The gauge field is incorporated in a real spinor superfield which, in the Wess-Zumino 
gauge, takes the form 

V a = t(A0) a + ^99 Va (86) 

where r) a is the supersymmetric partner of the gauge field (gaugino). 

The supersymmetric gauge invariant lagrangian for the matter fields which in terms 
of superfields is the highest component of the superfield 



VZ*VZ 



(87) 



with 

V a = D a - iV a (88) 
In terms of component fields the lagrangian reads: 

L = gl[D^z a D^z a + ill p^l + F a F a + 2i(r]^ a z a - W-qiP)] (89) 

where denotes the gauge covariant derivative with respect to the UsiX) field, and for 
convenience we have rescaled the fermion fields \l/ and the auxiliary field F by g%, as 
compared to the non-supersymmetric case, in order to facilitate our superfield formalism. 

Notice that ( |8~9"D contains a supersymmetric partner (gaugino) of the statistical gauge 
field Us(l)- This defines the N — 1 supersymmetric point of the model, in the sense that 
the gauge interaction Us(l) 'doubles' its degrees of freedom as a result of supersymmetry. 
The interactions of the gaugino rj wth the matter fermion (holon) \l/ and its partner, the z 
magnon (spinon), lead to an effective electric-charge violating interactions on the spatial 
planes, given that the Majorana gaugino rj is a real field, and as such cannot carry electric 
charge (which couples as a phase to a Dirac field). These terms can be interpreted as the 
removal or addition of electrons due to interlayer hopping ]T3 1 . 



It is important to notice that the constraint Q8T| ) admits a N = 1 supersymmetric 
formulation, in terms of the superfields Z a : 

jZZ a Z a = l (90) 

a=l 

which in components yields the constraint (]81f) as well as [13, 23]: 

Zj> a = (91) 

This can be solved by means of a 'colourless' fermion field X that satisfies (on account of 

= e a pZpX , X = EapZo^p (92) 

The auxiliary fields F a can be solved by means of their equations of motion and the con- 
straint (0), or, equivalently, in a path integral formalism by implementing the constraint 
via a Lagrange multiplier superfield: 

A(x, 6) = a(x) + 98{x) + -60\(x) . (93) 

In the second method, by eliminating F one obtains F a = —az a , for each a. If one uses 
the bosonic part of the super-constraint (|9~0|) , J2a=i \ z a\ 2 = 1, one then obtains: 

^« = E J 2 ( 94 ) 



We therefore observe that the supersymmetric extension of the CP 1 o model necessitates 
the presence of attractive Gross- Neveu type interactions among the Dirac fermions of each 
sublattice, in addition to the gauge interactions. An important point to notice is that the 



Gross-Neveu terms are of the type that would violate parity if dynamical generation of 
fermion mass occured as a result of these interactions. 

For completeness, we also note that the presence of Chern-Simons terms, which may 
appear in the effective action as a result of parity violation, does not add any complciation. 
As discussed in ref. f23|, the supersymmetrization of these terms leads to a mass for the 
gaugino, as expected from the fact that the Chern-Simons term is a topological gauge 
boson mass term J34|. The result (in components ) is: 



Qsupersymm 
^CS 



tuv p aud u a p + -7)7) 



(95) 



where 7) is a Majorana fermion, as we mentioned before, and k denotes the coefficient of 
the Chern-Simons term. 

Extension to N = 2 Supersymmetric CP 1 models 



The supersymmetric N = 2 CP M a- mo del was constructed in ref. |2B[], by dimensional 
reduction from a four- dimensional N = 1 supersymmetric lagrangian in a Minkowskian 
space time, which is super gauge and U(M) invariant: 
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d A xd 2 ed 2 6<&e v <$> 



d*xd 2 6d 2 6VDDV 



(96) 



in a standard four- dimensional superfield notation [35| with 6 a , 9 a complex spinors, and 
D a = gfn' ~ w^"^, D° = -g=- + i9 a cr^d '^j the vector superfield V is real, and contains 

the gauge bosons, A^, whilst the scalar superfield $ is chiral (D a § = 0). In component 
form the action (|96f) reads: 



5-( 4 ) = i j d 3 xdt 
i J d 3 xdt 



D^zD^z + itp ptp + FF + Xzz + ?puz + uipz 



+ 



1 

T 2 



-F 2 + 2iu Bu + 2d 2 



irrelevant 



(97) 



where D„ = cL — iAu, and the various component fields can be understood by the field 



content of the supersymmetry algebra PS] . The terms marked irrelevant are, by power 
counting, irrelevant operators in a renormalization groups sense in the low-energy regime 
of the dimensionally-reduced three-dimensional theory, and hence they can be safely ig- 
nored. The dimensional reduction in the Minkowski time t leads to a Euclidean three 
dimensional theory, which is precisely a N = 2 supersymmetric u-model [11]: 



S 



(3) 



DizD { z + itp ptp - A\zz + A Q ipip - FF + Xzz + ipuz + ui>z (98 



up to irrelevant operators of the form -p \F 2 ^ — ^(diA ) 2 — 2%u jdu — 2X 2 . Notice that 
the temporal component of the (dimensional reduced) four-dimensional gauge potential 
A Q plays the role of the Lagrange multiplier a field in an N = 1 formulation (c.f. 



The important differnce of the N = 2 formalism, however, is that now the gaugino field 
u is a Dirac spinor. 

A N = 2 supersymmetric version of the Chern-Simons terms also exists [23]. In 
component terms is given by 



qN=2 supersymm / t3 

D C5 — a l 



■CCKj 



1 



(99) 



The mixing Act is a feature of the N = 2 formalism and was absent in the N — 1 case 

(El)- 

In the context of our statistical spin=charge separating model, the presence of a Dirac 



gaugino allows, in contrast to the N — 1 case, for electric charge conservation |13j, given 
that now the gaugino being a Dirac spinor is allowed to carry electric charge. Thus the 
coupling terms to the matter fermions are now electrically neutral. This would imply 
suppression of intersublattice or interlayer hopping. 

The critical behaviour of the N = 2 model has been studied in detail in ref. . Due to 
the extra supersymmetry, the N = 2 case allows for an exact computation of the pertinent 



renormalization-group /5-functions [f^, [lT| , and hence the critical exponents of the model. 
This has to be contrasted with the situation in the N — 1 model, where such exact 
results are not available. Such exact results make the N = 2 CP M cr-model attractive for 
further studies along the lines of the fully non-perturbative approach to strongly-coupled 
supersymmetric gauge theories, advocated by Seiberg and Witten |L2|, [14]]. Then, by 
viewing the N — 1 case, of relevance to us in the context of doped antiferromegnets, 
as a supersymmetry boken descendant of the N = 2 model, one might obtain valuable 
non-perturbative information for the phase structure of the theory at the supersymmetric 
points of the parameters of the system. 
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